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Abstract. We investigate the effect of Fermi motion on the direct fragmentation of the J/xj) and 
X states employing a light-cone wave function. Consistent with such a wave function we set 
up the kinematics of a heavy quark fragmenting into a quarkonia such that the Fermi motion 
of the constituents split into longitudinal as well as transverse direction and thus calculate the 
fragmentation functions for these states. In the framework of our investigation, we estimate 
that the fragmentation probabilities of J/ip and X may increase at least up to 14 percent when 
including this degree of freedom. 
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1 Introduction 



Quarkonia production in high energy hadron colli- 
sions have shown to be an important issue both in 
QCD and in collider physics. Hadronic production of 
J I ip and T have received attention since the discovery 
of charm and beauty quarks. The production of these 
states have been studied in interesting models [1]. It 
is also found that there is an order of magnitude dis- 
crepancy between cross section data and theoretical 
calculations in the case of J/i[> at the Tevatron ener- 
gies. Then it became clear that the dominant mech- 
anism at high energy collisions are the color singlet 
[1,2] and the color octet [3] mechanisms and more 
over that at the Tevatron energies the color singlet 
contribution is too small [4,5]. In the color singlet 
mechanism it is assumed that the QQ pair is pro- 
duced in color singlet state at the production point 
whereas in the later color octet case the QQ pair is 
produced in the form of color octet. Later on emission 
of a soft gluon takes the pair into a color singlet state. 
To produce singlet or octet QQ there will of course 
be various channels out of which direct quark and 
gluon fragmentation are thought to be the dominant 



ones. In each case we need fragmentation functions 
to evaluate the probabilities and cross sections. 

The very accurate and successful calculation of 
the fragmentation functions have recently been car- 
ried out in the limit of heavy quark effective theory. 
It has been proved that the strong interaction per- 
turbation theory is suited to calculate the hard part 
of the heavy meson fragmentation in the above limit 
[6] and use a kind of wave function to introduce the 
soft smearing of the bound state. These are then con- 
voluted to give the fragmentation function. This con- 
volution is possible only in a certain QCD scale. The 
Altarelli-Parisi evolution equation [7] takes this func- 
tion into a required scale. In the case of fragmentation 
of heavy mesons, either a kind of delta function or the 
wave function at the origin of the bound state have 
been employed which produces a factor of meson de- 
cay constant, Jm, or the wave function at the origin 
i?(0)| and lets the constituents move almost together 
in parallel ignoring their virtual motion. However the 
problem of J/tJj cross section at Tevatron mentioned 
earlier and also its diffractive production [8] require 
a closer look at the process. 

The common procedure for introducing the bound 
state effects, is to use the wave function at the center 



2 



M.A. Gomshi Nobary and B. Javadi: Direct Fragmentation of Quarkonia.. 



of momentum frame in a harmonic oscillator model. 
However the kinematics and parameterization in such 
a case is not much clear. On the other hand the pro- 
cedure has successfully been studied in an infinite 
momentum frame. And that Brodsky-Huang-Lepage 
(BHL) [9] have proposed a boost free wave function 
for the heavy meson bound state which fits in the 
above calculational scheme. This is a kind of gaus- 
sian wave function with a parameter which adjusts 
its width is referred to as the confinement param- 
eter. The mass of the constituents and their trans- 
verse momentum as well as the longitudinal energy- 
momentum ratios appear in this function. 

In this work we first describe the BHL wave func- 
tion. Using it we calculate the fragmentation function 
for the J/tp and Y states in both cases i.e. when the 
Fermi motion is either ignored or not ignored. Change 
of the confinement parameter provides a nice com- 
parison of these two functions. It is revealed that the 
inclusion of this degree of freedom may have a siz- 
able effect on fragmentation probabilities and hence 
on their cross section measured in a particular ma- 
chine aimed to produce such states. 



2 Heavy meson wave function 

The light cone wave functions described in this sec- 
tion have satisfactorily been applied for light meson 
states such as ir, p and K [10]. motivated by these 
studies we will apply them to the case of J/tp and T 
states. 

The mesonic wave function can be obtained by 
solving the bound-state equation [11] 



where 



(M 2 ~H LC )\M) = Q, 



(1) 



where Hlc is the hamiltonian in light-cone quanti- 
zation and M is the mass of the meson \Ai). In an 
attempt to solve this equation one expands \M) in 
the complete set of Fock states 



1-^) = ^2 t l ; n(x l ,q Ti ,X l )\n;x i ,q Tl , A;}. (2) 

n,\i 

In the above relation ip n (xi, , Aj), (i = 1, 2, n) is 
the n-parton momentum-space amplitude defined on 
the free quark and gluon Fock basis in the light-cone 
formalism. The normalization condition is 



n n 

[dx] ee Y[ dx-iS 1 - Xjj , 



i=l 



i=l 



and 



n n 

[rf 2 q T ]^n rf2 ^ 16 ^ 2 [E^ 



1=1 



(4) 



(5) 



The sum is over all Fock states and helicities. 

It is not easy to determine the form of the wave 
functions from the above expansion. The prescrip- 
tion due to Brodsky-Huang-Lapage (BHL) [9] pro- 
vides a practical way by connecting the equal time 
wave function in the rest frame and the light-cone 
wave function, i.e. by equating the off-shell propaga- 
tor e = M 2 — (X)"=i 1i) m two frames: 



i=i 

n 

= M 2 -]r[(q^+m 2 )M] [L.C.] (6) 
i=i 

where the constraints ^™=i Pi — ^ an d Ym=i QTi = 
0, Yl^—i Xi = 1 apply for cm. and L.C. frames respec- 
tively. It is easily seen that for a two particle system 
with equal mass {p\ = p^) 



q 2 



4x(l - x) 



m 



(7) 



where p is the moment of the constituents in the cen- 
ter of momentum frame. Therefore the wave func- 
tions at two frames may be connected by 



A.mip 2 ) < > IpLC 



r q 2 ^ + m 2 



4x(l - x) 



(8) 



Now if we consider the wave function at the rest 
frame as [12] 



V>c.m.(p 2 ) = ^Mexp 



2(i 2 



(9) 



]T J [dxp 2 q T ]|^,q Tt ,A 4 )| 2 = 1, (3) 



where (3 is the confinment parameter, then using the 
connection (8) the IMF wave function in the case of 
equal mass constituents takes the following form 
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Note that in this formalism = — qr 2 , therefore 
gij. = q^ 2 = q\- It is argued that this form holds 
in the case of heavy meson states with equal and 
unequal constituent quark and antiquark masses, i.e. 



ip M {xi,k T ) = A M exp< 



'IT 



L 2T 



8/3 2 L X\ X2 



.(11) 



q| and 



- q|n. Tfi\ and 



TO2 ar e quark or anti-quark masses chosen so that 
mi > rri2- 

It is interesting to note that since the spin of the 
heavy quark decouples from the gluon field, the ex- 
ited states has the same wave functions as the un- 
excited ones. The behavior of ip M is shown in Fig. 
1 for the J/tp state with three different confinement 
parameter (3 which are employed to extract the effect 
of Fermi motion in this work. Naturally the shape of 
the wave function is very sensitive to the values of j3. 
x\ and xi both go from to 1 with the constraint of 
x\ + X2 = 1. qr should cover the necessary range to 
include the effective spreading of the wave function. 



13 = 0.2 




[3 = 0.4 




(3 = 0. 




3 Kinematics 

Typical lowest order Feynman diagram for fragmen- 
tation of a c(b) quark into J/ip(T) is shown in Fig. 2. 
To evaluate such a diagram consistent with the light 
cone wave function for the bound state, we need to 
choose an infinite momentum frame in which the ini- 
tial state heavy quark possesses transverse momen- 
tum of p't- This transverse momentum is to be car- 
ried by the final state c(6) which produces a jet. The 
constituents of the meson after creation, move along 
the z axis with transverse momentum of pt = — 
(pt = &t = <zt)- Therefore the particles' four mo- 
menta are considered as 

Vn = PL, Qt ] , k^= [fc , k L , -q T ] , 
p'^= [p' , p' L , k T ], k'^[k' , k' L ,k T }. (12) 

Note that all particles are moving in the forward di- 
rection and qr's are chosen in a consistent manner 



Fig. 1. The behavior of the BHL meson wave function 
(11) squared for the J/ip state with different values of 
the confinement parameter, x is the longitudinal energy 
momentum ratio and q is the magnitude of the trans- 
verse momentum of the constituents. The scale for ip 2 is 
arbitrary. 



with the meson wave function. Moreover if we desig- 
nate the meson four momentum by p, then we may 
write p M = Xip^ and fc M = X2p^ with the condition 
that x\ + X2 = 1- Therefore, we parameterize the 
energies of the particles as 

p = xizp' Q , k Q = x 2 zp' , k' = (l- z)p' , p' Q = pf^i) 

where we have used the definition of the fragmenta- 
tion parameter as z = -Ehadron/^-quark in an infinite 
momentum frame consistent with our kinematics. 
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Q(p') 




(p a + k Q + k' Q — p' Q ) is the energy denominator. In the 
case of Fig. 2, r has the following form 



r = u(p)j^u(p')u(k')j tl v(k) . 



(15) 



The initial spin average and final state spin sum of 
the squared r takes the following form 



1 



( P .k')( P '.k) + (p.k){p'.k') 



+m%(p.p') - ml(k.k') - 2m\m 2 2 .(16) 

Regarding the kinematics introduced earlier, the 
dot products in the above read 



Fig. 2. The lowest order Feynman diagram for fragmen- 
tation of a heavy quark into a quarkonium state. The four 
momenta are labeled. 



4 The process of fragmentation 

The very complex process of a heavy quark fragment- 
ing into a heavy mesonic bound state is understood 
and becomes also calculable if we take into account 
the following 

• The hard part of the process of a heavy quark 
splitting into a heavy meson is calculable in the per- 
turbative QCD in fragmentation scale [13]. The wave 
function of the bound state is convoluted to the hard 
scattering amplitude in an appropriate manner. 

• To absorb the soft behavior of the bound state 
into hard scattering amplitude a suitable wave func- 
tion is employed (our choice (11)). 

The hard scattering amplitude for the creation of 
the bound state may be put into the following form 
[14] 



T H = 



2^2p p' k k' 

r 



{k + k'Y{ Po + k + k' - P ' y 



(14) 



2p.k' = -(m\ + g|) + p— (m\ + k 2 T ) 

X\Z 1 — z 

-2q T .k T , (17) 



2p '.k = —(ml + q T ) + (m? + k 2 T ) - 2q T .kft8) 

XoZ 1 



2p.k = ^-(m% + 4) + — (m\ + 4) + 2q%, (19) 

X2 i 



2p'.k' = l —^(m\ + k%) + — !— (ml + kl) - 2fc^0) 
1 1 — z 



2k. k' = - — -(ml + qt) + ^—(m 2 2 + ki) 

X2 Z 1 — Z 



-2q T .k r , 



(21) 



2p'.p =—(mj+ 4) + ^(mj + k 2 T ) - 2q T .k£22) 

X\Z 1 

Here qy-kp = qr-kr cos where 9 is the angle be- 
tween qr and kp. 

The fragmentation function of a heavy quark frag- 
mentation into a heavy meson with wave function ip 
is obtained from [14] 



Here a s = ,g 2 /47r is the strong interaction coupling 
constant, Cf is the color factor and r indicates that 
part of the amplitude which embeds spinors and gamma 
matrices, (k + k') 2 is due to the gluon propagator and 



D{z,ii ) = ^J2 f d3 pd 3 kd 3 k' 

2 s J 

x\T H \ 2 \^ M \ 2 S {3) (p + k + k' - p').(23) 
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Here Th is the hard scattering amplitude given by 
(14) and ipM is the meson wave function (11). Note 
that average over initial quark spin and sum over final 
state particles spin are assumed. We replace (14) into 
(23) to obtain 



D(z,fj, ) = 



{4ira s m 1 m 2 MC F ) 2 f d 3 pd 3 kd 3 k' 



PoP k k 

^Tr\^ M \ 2 s^( P + k + k' - v 'l 



f{zf{ Po + k + k' - P ' Y 



where f(z) — (k + k') 2 . Now we perform the phase 
space integrations. First of all note that 



/ 



d 3 P 5 {3 \p + k + k' - P ') 
p'Jpo + k + k' Q - p' a ) 2 ~ g{z) 2 ' 



P'o 



(25) 



where g(z) = (p + k + fc') 2 When the Fermi motion 
is ignored, \^m\ 2 is integrated with d 3 k and gives 



/ 



d 3 £# M | 2 = 



ZPo 

16tt 3 ' 



(26) 



where we have used the normalization condition (3) 
on the wave function. Also note that in this case we 
have 



J d 3 k' = p' Q J d 2 k' T = 2np' J k' T dk' T . (27) 

Instead of doing this integration, we simply replace 
k' T by its average value. Thus we obtain 



note that we have distinguished the quark and anti- 
quark masses such that (28) and (29) may be used in 
the case of other meson states such as D, B and B c . 

The explicit form of the fragmentation functions 
(28) and (29) in general are rather lengthy. However 
in the case of quarkonia they cast into rather brief 
forms. They appear in the appendix. 



5 Results 

We have shown that the wave function due to BHL 
may be employed to represent the soft behavior of 
the heavy meson bound state. In infinite momentum 
frame we have set up a kinematics which embed the 
Fermi motion of the constituents in both longitudinal 
and transverse direction through this wave function 
whose width is controlled by the confinement param- 
eter (3. Comparison of (28) and (29) demonstrates 
the effect in an illustrative manner. This compari- 
son is shown in Figure 3 in the case of the J/tp and 
Y states. The two functions coincide at sufficiently 
small values of (3. As we increase (3, the peak due to 
(28) rises and (z) moves to higher values. This be- 
havior is consistent with the physics of the effect. In 
fact inclusion of Fermi motion in the fragmentation 
of a particular state means fragmentation of a state 
with higher energy and momentum which naturally 
ends up with higher values of both the fragmenta- 
tion probability and average fragmentation parame- 
ter. Note the ranges of (3 selected in the two diagrams 
in this figure. They are the same as those used in Fig- 
ure 1. To find the size of such changes, we employ the 
fragmentation probability defined by 



D(z,Mo,/? = 0) = 



mi {a 8 MCF) 2 (k' T 2 ) 1/2 

r 

X z{l-z)f{zfg{zf 



(28) 



where /' represents the square bracket in (16). When 
the Fermi motion is included, along the same line of 
calculation we obtain 



D(z,fj, ,P) = 



m 1 {^a s MC F ) 2 {k' T 2 ) 1 / 2 



dq T dx 2 qTf'\^M\ 2 
x 1 x 2 z(l - z)f(z) 2 g(z) 2 



.(29) 



Note that the above scheme of the phase space in- 
tegration allows one to compare (28) and (29) sim- 
ply by varying the confinement parameter (3. Further, 



P(0)= J D(z,^ 0l P)dz, (30) 
and the average fragmentation parameter 



(*)(?) = 



J zD(z,Ho,P)dz 
J D(z,n ,f3)dz ' 



(31) 



We plot P{(3) I P(0) versus (3 in Figure 4. Such evalua- 
tions reveal that the Fermi motion of the constituents 
of J/ip and Y states improve the fragmentation prob- 
abilities of these states at least up to 14 percent. We 
have also plotted {z)((3) against (3 in Figure 5. 

This investigation shows that the Fermi motion 
within the heavy meson states such as J/4> and Y 
is indeed important and could not be ignored. This 
of course will be reflected on the cross section and 
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Fig. 3. The affect of J/tp and T wave function on their 
fragmentation. While the dashed curve represent (28) , the 
solid curves show the behavior of (29) with increasing val- 
ues of (5. The two curves coincide exactly for sufficiently 
small values of confinement parameter. 
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Fig. 4. Rate at which the fragmentation probabilities for 
J/tp and X increase with increasing confinement param- 
eter. Note the ranges of f3 selected for the two cases. 

event rate of these states at the colliders where they 
are expected to be produced. 

Concerning the J/tp production our results could 
be accounted for as a modification to the color singlet 
mechanism where contribution due to various chan- 
nels specially the color octet mechanism are relevant. 
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Appendix 

The explicit form of (28) and (29) for quarkonia cast 
into the following 

*lC%{k' T 2 )^ 
WmJ 7 

x {z(l - zf [£V + 2£z 2 (4 -Az + hz 2 ) 
+ (16 - 32z + 24z 2 - 8z 3 + 9z 4 )] j, 

n*a* s Cl(k' T 2 )V* 
2m 

f dqdxqx 2 z(l - z) 2 \i\j m \ 2 

X J g w ' 

where T , Q and H read as 

T = [fz 4 -(z- 2) 2 (3z - 4) + £z 2 (8 -7z + z 2 )] \ 



Q = { [77(1 - z) 2 + ix 2 z 2 + (1 - (1 - x)zf\ 

x [r](-l + z) + £(-1 + x)xz 2 - 1 + (1 - x + x 2 )z)] y, 

and 

H = |l - 4(1 - x)z + 2(4 - 10a; + 7a; 2 )z 2 + 4(-l + Ax 

-5a; 2 + x 3 )z 3 + (1 - 4x + 8a; 2 - 4a; 3 + x 4 )z 4 + r]£z 2 
x [l - 2x + z 2 + x 2 (2 -2z + z 2 )] + 77(2 + (4a; - 6)z 
+ (9 - 8a; + 2x 2 )z 2 - 2(2 - x + x 2 )z 3 + (1 + x 2 )z 4 ) 
+£z 2 (l + 2x 3 (2 - 3z)z + z 2 + 2x 4 z 2 + x(-2 + 2z 
-Az 2 ) + a; 2 (2 - 8z + 9z 2 )) + r) 2 (l - z) 2 

+e 2 (l-.x)W}. 



Here £ = (k' T )/m 2 and 77 = q 2 /m 2 . Am is calcu- 
lated using the normalisation condition of the wave 
function, e.g. equatrion (3). 



D(z,Ho,0 = O) = 



and 

D(z,/jl ,/3) = 



